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Abstract. Young galactic X-ray point sources (XPs) closely trace the ongoing star formation in galaxies. From measured 
XP number counts we extract the collective 2-10 keV luminosity of young XPs, L^'^^ , which we use to gauge the current 
star-formation rate (SFR) in galaxies. We find that, for a sample of local star-forming galaxies (i.e., normal spirals and mild 
starbursts), LJ^^ correlates linearly with the SFR over three decades in luminosity. A separate, high-SFR sample of starburst 
ULIRGs can be used to check the calibration of the relation. Using their (presumably SF-related) total 2-10 keV luminosities 
we find that these sources satisfy the SFR-LJ'^^ relation, as defined by the weaker sample, and extend it to span ~5 decades 
in luminosity. The SFR-LJ'^^ relation is likely to hold also for distant (z~l) Hubble Deep Field North galaxies, especially so 
if these high-SFR objects are similar to the (more nearby) ULIRGs. It is argued that the SFR-LJ^^ relation provides the most 
adequate X-ray estimator of instantaneous SFR by the phenomena characterizing massive stars from their birth (FIR emission 
from placental dust clouds) through their death as compact remnants (emitting X-rays by accreting from a close donor). For 
local, low/intermediate-SFR galaxies, the simultaneous existence of a correlation of the instantaneous SFR with the total 2-10 
keV luminosity, Lx, which traces the SFR integrated over the last ~10'"* yr, suggests that during such epoch the SF in these 
galaxies has been proceeding at a relatively constant rate. 

Key words, galaxies: X-rays - galaxies: spiral - galaxies: starburst - galaxies: stellar content - galaxies: evolution - infrared: 
galaxies - radio continuum: galaxies - stars: formation - X-rays: binaries - X-rays: galaxies 



i 1. Introduction 

Star formation (SF) leads to X-ray emission on various spatial 
and temporal scales, including emission from O stars. X-ray 
• binaries, supernovae (SN) and their remnants (SNRs), galactic- 
scale emission from diffuse hot gas, and Compton scattering of 
' FIR & CMB photons by SN-accelerated electrons. Integrated 
' spectra of star-forming galaxies (SFGs) are expected to show 
all these components, as well as emission from an active nu- 
cleus (Rephaeli et al. 1991, 1995). Based on lower-resolution 
X-ray data (mainly from Einstein, ASCA, BeppoSAX, and 
RXTE), most SFGs have remarkably similar integrated spec- 
tra that include a soft (single- or multiple-temperature) sub- 
keV thermal component which dominates at energies 2 keV, 
plus a hard power law (PL) which dominates at e<; 2 keV (e.g., 
Dahlemetal. 1998). 

By quantitatively assessing the spectral components of 
the various X-ray emission mechanisms in SFGs, Persic & 
Rephaeli (2002) concluded that the 2-10 keV emission is dom- 
inated by compact X-ray binaries, specifically HMXBs if the 
star-formation rate (SFR) is very high. 

Stellar-related X-ray emission can be used as an indicator 
of SFR (David et al. 1992; Bauer et al. 2002; Grimm et al. 2003; 
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Franceschini et al. 2003; Ranalli et al. 2003; Gilfanov et al. 
2004a). The basic notion is that the ongoing SFR can be mea- 
sured based on stellar end-products which are both sufficiently 
X-ray-bright for their collective emission to be unambiguously 
identified, and sufficiently short-lived so that they trace the 'in- 
stantaneous' SFR. Of the three main types of Galactic stellar 
endproduct X-ray sources (LMXBs, SNRs, and HMXBs), the 
latter two provide a suitable combination of short delay be- 
tween star formation and onset of X-ray emission and signif- 
icant X-ray brightness. 

Persic et al. (2004a) examined ways in which the 2-10 keV 
luminosity (hereafter: Lx) can be used as an estimator of ongo- 
ing galactic SFR. They concluded that the collective 2-10 keV 
emission from HMXBs, L™-^^, was a reliable SFR estima- 
tor (supporting independent suggestions by, e.g., Grimm et al. 
2003). Given the diverse sources of X-ray emission in SFGs, 
Persic et al. (2004a) argued that the level of HMXB emission 
could be evaluated by modelling SFG spectra with Persic & 
Rephaeli's (2002) template - in which the HMXB component 
was represented as a r=1.2 power law. Upon analyzing suit- 
able ASCA, BeppoSAX, and Newton/XMM spectra, Persic et al. 
(2004a) suggested that L™™-- 0-2 L.^ in moderately star- 
forming galaxies and L^^^^^^Lx in intensely star-forming 
ones. 



2 Persic & Rephaeli: Stellar birth-rates from XPs 

Table 1. Data la: local normal and starburst galaxies (SFG sample). 
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Distances (in Mpc) are from Tully (1988) if cz < 3000 km s ^, or consistently assume Ho = 75 km s ^ Mpc ^ otherwise. 

Integrated 2-10 keV luminosities (in erg s^^; given in log form). For references see Persic et al. (2004a). For NGC 3077 and NGC 4214 see 
below. For NGC 4945 see Guainazzi et al. (2000). 

Cumulative XP luminosities in the 2-10 keV band (in erg s~^; in log form). The quoted values result from converting the values measured 
in several different bands [0.3-10 keV for NGC 2403 (Schlegel & Pannuti 2003); 0.5-5 keV for NGC 6946 (Holt et al. 2003); 0.2-12 keV for 
IC 342 (Kong 2003); 0.5-8 keV for Arp 299 (Zezas et al. 2003); 2-8 keV for NGC 891 (Temple et al. 2005); 0.3-8 keV for NGC 4214 (Hartwell 
et al. 2004) and for all other objects (Colbert et al. 2004)] into the 2-10 keV band through an assumed photon spectrum cj) oc e~" [with a = 2 
in Kong (2003), 1.8 in Hartwell et al. (2004) and Temple et al. (2005), and 1.7 elsewhere]. 

In a few cases (NGC 4038/9, NGC 4449, NGC 5457, NGC 6946) it turns out that logLxp > logL^ by 0.10 < 6logL< 0.30: assuming such 
unphysical situation to descend from non-optimal fits being performed on most XP spectra (except for those few XPs in each galaxy which had 
enough counts to allow detailed fits to be performed), in our computations we set Lxp=ia:. 

For NGC 4214 the XPLF slope is taken slightly steeper than quoted in Hartwell et al. (2004), to account for incompleteness that admittedly 
affects the XPLF shape below logL = 36.6. 

For NGC 2403 the reported value of ixp is based on the mean flux from its XPs (Schlegel & Pannuti 2003), excluding the 4 brightest sources 
and 12 presumed interlopers. The reported 7 and 77 refer to the same set of sources represented by Lxp. See discussion in Section 4.1. 
For NGC 2146 the 2-10 keV XP data are taken from Inui et al. (2005). 

For NGC 3077 the luminosities of its 6 XPs and of the diffuse thermal plasma, measured in the 0.3-6 keV band (Ott et al. 2003), have been 
converted to 2-10 keV luminosities, using the individual published spectral models of Ott et al. (2003), and combined to give L^. The reported 
Lxp refers to the estimated young-XP luminosity. No published XPLF is available. See discussion in Section 4.1. 

^"^^ Minimum and maximum XP luminosities (in erg s~ ^ ; given in log form), between which the XPLF index, 7, was derived. For spectral band 
information, see previous point. 

Index of the cumulative XPLF, of the form N{> L)ocL ''. For references, see point (c). 
^■^^ Blue apparent magnitudes, corrected to face-on and for Galactic absorption, from RC3 (de Vaucouleurs et al. 1991). 

IRAS flux densities at 12/im, 25Atm, 60/im, and lOO^m (in Jy), from Sanders et al. (2003; except for NGC 4449: Hunter et al. 1986). The 
FIR flux is accordingly defined (Helou et al. 1985) as /fir = 1.26 x 10^^^(2.58/60 + /ioo)ergs"^ cm"\ 

1.4 GHz flux density (in Jy). Data are from Kiihr et al. 1981 (NGC 253), White & Becker 1992 (NGC 628, NGC 891, IC 342, NGC 2403, 
NGC 3034, NGC 3079, NGC 3628, Arp 299, NGC 4631, NGC 5457, NGC 6946), Condon et al. 1998 (NGC 1569, NGC 2146), Heeschen & 
Wade 1964 (NGC 4449), Wright & Otrupcek 1990 (NGC 4945, NGC 5236). No entry means that no data, or only upper limits to /1.4, were 
found in the literature. 

Notes, (i) Compared with Persic et al. (2004a), several objects could not find their way into the present analysis. These systems and the reason 
for flieir exclusion are: NGC 55, NGC 891, NGC 1808, NGC 2276, NGC 5782, NGC 2903, NGC 3256['l , NGC 3310''1 , NGC 3367, NGC 3556, 
NGC 4654, NGC 4666''' , NGC 7552: no puWished XPLF available ([/] : observed with Chandra and/or XMM-Newton (Lira et al. 2002; Jenkins 
et al. 2004; Persic et al. 2004b). NGC 5253: no integrated 2-10 keV luminosity available (see Summers et al. 2004). NGC 7469, NGC 7679: 
AGN-dominated 2-10 keV emission, (ii) Other commonly used names for some of the objects are: M 74 for NGC 628, M 82 for NGC 3034, 
NGC 3690 + IC 694 for Arp 299, Antennae for NGC 4038/9, M 83 for NGC 5236, and M 101 for NGC 5457. 
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Table 2. Data lb: local normal and starburst galaxies (SFG 
sample). 
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Derived 2-10 keV luminosity fraction of young XPs. 
'''^ Cirrus-uncorrected FIR luminosities (in erg s^^), in log form. 

Cirrus-corrected FIR luminosities (in erg s^^), in log form. 

Star-formation rates, in A4q yr^^, uncorrected for cirrus emis- 
sion. They are derived from IR luminosities using eq.(3) and setting 
Lir=1.65xLfir. 

Star-formation rates, in M© yr"'^, corrected for cirrus emission. 
They are derived as described in point (d). 

Note. For objects located at low Galactic latitudes (\b\ < 15 degrees) 
the corrections for foreground Galactic absorption tend to be large and 
uncertain. An extreme case is IC 342 (&=10.58 deg), for which the B- 
band absorption is estimated to be as large as 3.360 mag (Burstein & 
Heiles 1982) or 2.407 mag (Schlegcl et al. 1998). As the adopted sta- 
tistical correction for cirrus emission turns out to be unphysical in this 
case, we will drop IC 342 from further analyses. None of the results 
reported in this paper will depend on the exclusion of this object. 



However, Persic et al.'s (2004a) approach was potentially 
limited. The template SFG spectrum of Persic & Rephaeli 
(2002), on which their spectral analysis was based, exphcitly 
assumed a Galactic stellar population: in principle that fea- 
ture could bias Persic et al.'s (2004a) procedure because any 
SF-tracing population of X-ray point sources (XPs), not repre- 
sented in our own Galaxy but present elsewhere (and spectrally 
different from Galactic HMXBs), would be missed. This ac- 
tually turns out to be the case for the Ultra-Luminous X-ray 
sources (ULXs) which dominate the XP luminosity (Lxp) of 
most SFGs. Furthermore, few integrated spectra were of suffi- 
ciently good quality to be fitted with the multi-component tem- 
plate of Persic & Rephaeli (2002). 

It is important then to further scrutinize the role and ef- 
fectiveness of the collective 2-10 keV emission of young XPs, 
L^^^, in measuring the ongoing SFR in galaxies - as well as 



ways to evaluate it. In this paper we propose to evaluate 
by modelling the luminosity functions of galactic XPs (here- 
after: XPLFs) as linear combinations of the 'universal' young- 
and old-XPLFs. To demostrate the effectiveness of the pro- 
posed technique, we will apply it to a sample of galaxies with 
available Chandra XPLFs. For the same sample, we will then 
check how the derived values of L^^^ correlate with the on- 
going SFR as deduced from the galaxies' (thermal) FIR and 
(non-thermal) radio emission. 

The plan of this paper is as follows. In section 2 we review 
the properties of XPLFs in galaxies. The two local and nearby 
samples of star-forming galaxies are described in section 3. In 
section 4 we review the FIR- and radio-based SFR indicators. 
In section 5 we apply our proposed young/old-XP decomposi- 
tion technique to a sample of SFGs and evaluate L^.^^, which 
are then contrasted with FIR/radio-based SFRs: the results are 
discussed in section 6. In section 7 we extend our investigation 
to high-SFR galaxies, and discuss SF in the nearby universe. In 
section 8 we discuss X-ray SFR indicators at high redshift. X- 
ray and radio SFR indicators are compared in section 9. Section 
10 closes with a short summary of our main results. Throughout 
this paper we assume fio = 1, Ho=75 km s~^Mpc~^. (No re- 
sult in this paper is substantially affected by choosing this one 
particular cosmology.) 

2. X-ray point sources in gaiaxies 

Observations of nearby galaxies led to the detection of XPs 
and their luminosity functions (XPLFs) down to limiting lu- 
minosities of ~10^^ erg s^^ (see Fabbiano 2005 for a re- 
cent comprehensive review.) Such XPs include SNRs, close 
binary systems in which the accreting object can be either 
a pulsar {L^ < 2 y. 10^^ erg s^^) or a black hole (BH) 
(2 X < L^/{cYg s-i) < 1039), and ULXs {L^>IQ^'^ erg 
s~^). The quoted limits correspond to Eddington luminosities 
for spherical accretion onto ~2 Mq and ~8 Mq BHs, respec- 
tively, which is the approximate range of BH masses that are 
thought to be therein attainable via ordinary stellar evolution. 

Young SNRs, which result from the explosion of short- 
lived ^ 5Mq progenitor stars and are X-ray bright for ^lO"* 
yr, trace current SF. X-ray binaries, in which X-ray emission 
results from mass accretion onto a compact stellar remnant 
(NS or BH) from a main-sequence donor, are of the high- 
mass (HMXB) type or the low-mass (LMXB) type according to 
whether the donor mass is 8 Mq or ^Mq (see Persic 
& Rephaeh 2002 and references therein). HMXBs and LMXBs 
therefore trace, respectively, the current or average past SFR. 
The accreting object in ULXs is presumed to be a stellar (super- 
stellar) mass BH, accreting at super- (sub-) Eddington rates. 
Association with intense SF activity suggests that most ULXs 
are of the HMXB type (e.g., Zezas et al. 2002) - their opti- 
cal counterparts having sometimes been positively identified as 
O stars (e.g., Liu et al. 2002). When observed in E/SO galax- 
ies (only occasionally, and Umited to Lx~ 2 x lO^^ gj-g s"^), 
ULXs are likely of the LMXB type (Sarazin et al. 2000; Kim & 
Fabbiano 2004). In both cases, ULXs appear to extrapolate X- 
ray binaries to higher accretor masses (e.g., Swartz et al. 2004). 
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The association of XPs with their host environment proves 
useful to measure XPLFs for uniform source population. 
Functions measured in starburst environments yield the young- 
XPLF, and XPLFs measured in 'sterile' environments (e.g., 
in elliptical galaxies) reproduce the old-XPLF. As measured 
in systems with high SFR (e.g., in NGC 4038/9: Zezas & 
Fabbiano 2002; see discussion in section 4.3), the differential 
young-XPLFs can be described by 



dL 



- = Uyfl L 



-/3y 



1.5 



(1) 



with the cumulative counts Ny{>L)(xL^'-J^^^\ In the follow- 
ing we shall take eq.(l) to represent the 'universal' young- 
XPLF. The statistically superposed functions from a sample 
of E/SO galaxies result in a completeness-corrected differential 
old-XPLF which is given by 



diVv 
dL 



= »^v,o X 



Li < i < -Lbr 

LhT ^ L < L2 



(2) 



with /3v,i~1.8 and /3v,2~2.8 the faint-end slope and bright-end 
slope, Xi^SxlQ-^^ erg s~^ and _L2^2xlO^^ erg s~^ the lim- 
iting luminosities, and Lbr'^SxlO'^^ erg s~^ the break lumi- 
nosity (Kim & Fabbiano 2004; see also Gilfanov 2004). The 
corresponding cumulative function has a low-L slope ~1 and 
a high-L slope 1 .8. In the following we shall take eq.(2) to rep- 
resent the 'universal' old-XPLF (see discussion in section 4.3). 
The break may highlight a change in the nature of the XP pop- 
ulation: as its value approaches the Eddington limit for an ac- 
creting NS, the break may signal the NS to BH transition in the 
LMXB population. 

Galaxies with mild ongoing SF, like our own Galaxy, have 
both young and old XPs. A direct separation of young and 
old XPs has been possible only in few cases (e.g., in M 81: 
Tennant et al. 2001). Usually, measured XPLFs reflect galaxy- 
integrated counts and hence result from a combination of 
young-XPLF and old-XPLF, whose relative normalization de- 
pends on the current to average-past SFR (see Grimm et al. 
2003; Gilfanov 2004). Indeed, based on 1441 XPs detected in 
32 nearby galaxies with different levels of SF activity, Colbert 
et al. (2004) concluded that XPs are linked to both the old 
and young stellar populations or, equivalently, to both the past 
and present SF activity (confirming the earlier suggestions of 
Fabbiano & Trinchieri 1985 and David et al. 1992). To esti- 
mate the current SFR, young XPs have to be culled out from 
the total population. 

3. Galaxy samples 

Sample 1 (see Tables 1,2) consists of local SFGs with available 
IRAS FIR fluxes, ASCAJBeppoSAX/RossiXTE 2-10 keV fluxes, 
and Chandra XPLFs. It spans ~3 decades in SFR (estimated 
from LpiR, see section 4), from levels typical of quiescent iso- 
lated spirals (^1 M0yr^^) all the way up to strong merging 
starbursts (~5OM0yr~^). As such, sample 1 is fairly repre- 
sentative of the range of SF activity in the local Universe. 

Sample 2 (see Table 3) comprises a set of nearby Ultra- 
Luminous Infra-Red Galaxies (ULIRGs) with available IRAS 
flux densities and XMM 2-10 keV fluxes. As suggested by 
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Fig. 1. IRAS color-color diagram of the galaxies of samples 1 
(fiUed dots) and 2 (empty squares). Areas (a) through (d) denote 
successive phases of starburst evolution according to Vega et al. 
(2005; respectively: early-, peak-, evolved-, post-starburst). 



their LpiR, these galaxies are sites of very intense star forma- 
tion (SFR>100 Af0yr^^). Their 2-10 keV spectra show no ev- 
idence of AGN emission and are reminiscent, in their shapes, 
of Galactic HMXBs: this suggests that the entire of these 
galaxies may be related to ongoing SF (see Franceschini et al. 
2003). This set of starburst ULIRGs probes the peak of SF ac- 
tivity in the nearby Universe. 

It is instructive to check the evolutionary stage of the SFGs 
represented in samples 1 and 2. As the spectral region ~8-120 
fxm (sampled by the IRAS broad-band filters), which is charac- 
terized by continuum emission from hot dust, is very strongly 
affected by heating processes associated with SF, the IRAS 
color-color plot can be interpreted in terms of evolution of the 
SF activity. Based on advanced spectro-photometric modeling, 
Vega et al. (2005) suggested that the evolution of a starburst 
can be described as a sequence of four main phases (see Fig. 1): 
(fl) an early-starburst phase, when the newly formed stars are 
still deep inside their placental clouds, and the escaping radia- 
tion field has not reached its peak emission; ffoj a peak- starburst 
phase, when most massive stars are produced and are still em- 
bedded in their placental clouds, the SED is dominated by hot 
dust emission, and the starburst reaches its hottest colors; fcj an 
evolved-starburst phase, when the current SFR has decreased 
dramatically, the young hot stars have emerged from their pro- 
genitor clouds, and the cirrus emission is important; and (d) & 
post-starburst phase, when the current SF is mainly due to the 
quiescent disk and the colors are evolving towards those of nor- 
mal spirals. According to this scheme, the location of our sam- 
ple objects in Fig.l suggests that the ULIRGs are dominated 
by a starburst in its peak, while most local SFGs represent later 
phases, from evolved- through post-starburst to quiescent. 



4. FIR, radio galactic-SFR indicators 

A significant fraction of the bolometric luminosity of an ac- 
tively star-forming galaxy is absorbed by interstellar dust and 
re-emitted in the FIR band. As the absorption cross-section of 
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dust is strongly peaked in the UV which inherently traces mas- 
sive SF, the FIR emission can be a sensitive tracer of the current 
SFR. As discussed by Kennicutt (1998a,b), there probably is 
no single calibration that applies to all galaxy types. However, 
based on the main characteristics of our composite sample and 
guided by the general principle that the FIR emission should 
provide an excellent measure of the SFR in dusty starbursts, 
we shall adopt the most appropriate FIR SFR indicator for our 
current purposes. 

A wide range of conversion relations between SFR and 
LpiR are found in the literature, based on either a starburst 
model or an observational analysis (e.g.. Hunter et al. 1986; 
Meurer et al. 1997; Kennicutt 1998a,b). Using the continuous- 
burst model of 10-100 Myr duration of Leitherer & Heckman 
(1995) and the Salpeter (1955) stellar initial mass function 
(IMF) with mass limits 0.1-100 Mq, Kennicutt (1998a,b) de- 
rived a conversion relation appropriate for starbursts: 

SFR = — — r^^? r M0yr-i (3) 

2.2 X lO^Sergs-i 

with ~30% uncertainty. Here Lir refers to the full 8-1000 /im 
band: for starbursts with typical dust temperatures and emissiv- 
ities, however, most of the emission falls in the FIR (~40-120 
/xm) band, LFiR/iiR"^0.6 (Helou et al. 1988). Strictly speak- 
ing, the relation in eq.(3) applies only to young starbursts em- 
bedded in optically thick dust clouds. In more quiescent SFGs 
the assumptions underlying the relation in eq.(3) are not veri- 
fied. Among these, the dust optical depth is lower, and a colder 
"cirrus" component, originating from the heating of the ISM by 
the galactic UV emission that is powered mostly by old stars, 
will contribute to the total FIR emission. 

In this paper we adopt Kennicutt's conversion in eq.(3) as a 
FIR-based SFR indicator. The insight on the evolutionary status 
of the galaxies in samples 1 and 2, discussed in section 3, sug- 
gests some more accurate way of applying eq.(3). As suggested 
by their FIR colors (see Fig.l), all the ULIRGs of sample 2 
and some of the SFGs of sample 1 are young, dusty, optically- 
thick starbursts that presumably meet the assumptions under- 
lying eq.(3) and hence can be straightforwardly treated with 
it and their SFR can then be directly estimated from the ob- 
served LiR. For the remaining, milder SFGs of sample 1, which 
have evolved past the peak starburst phase, some of the as- 
sumption underlying eq.(3) are not valid, so in principle their 
SFR can not be directly estimated from the observed Lir. In 
particular, the FIR emission of these galaxies should be cor- 
rected for cirrus emission before being used as a SFR indica- 
tor. Since a detailed spectro-photometric modeling (e.g., Vega 
et al. 2005) of our sample galaxies is beyond our immediate 
scope, in this paper we follow David et al. (1992) in adopt- 
ing Devereux & Eales's (1989) statistical cirrus correction. The 
basic assumption is that the strong empirical FIR-radio corre- 
lation holding for SFGs (Helou et al. 1985), interpreted as a 
consequence of massive SF (see below), is blurred by the SF- 
unrelated cirrus FIR emission, -LpfR, most notably so at low lu- 
minosities where in fact some nonlinearity occurs (see Condon 
1992; Bell 2003). Setting L^I^ozLb (the blue luminosity be- 
ing a proxy for the galactic stellar content), then the SF-related 
FIR component is then Z/|fj^=iFiR-a; Lb : the FIR-radio cor- 



relation is linearized and optimized if a; = 0.14. The FIR lu- 
minosities in sample 1, although corrected according to this 
recipe for SFR-computing purposes, for simplicity will still be 
called Lfir (see Table 2). Finally, the FIR (~40-120/xm) lu- 
minosities need a correction for the wider bandwidth (~8-1000 
yum) required by Kennicutt's formula. Although such correc- 
tion clearly depends on the detailed spectral energy distribution 
of each object, we here adopt a statistical correction and set 



^IR 



= fL 



FIR 



/ = 1.65 



(4) 



to be used in eq.(3). This bandwidth correction is, strictly 
speaking, valid for starbursts with /go //loo ~ 1 ™d dust emis- 
sivity index equal to (Helou et al. 1988) and hence it may ap- 
ply only to the ULIRGs of sample 2 and to some strong SFGs of 
sample 1 ; however we assume that, after removal of the cirrus 
component, it is sensible also for the remaining, more mildly 
star-forming objects of Tables 1 and 2. In substantial agreement 
with our choice, Hopkins et al. (2003), Bell (2003), Kewley et 
al. (2002), and Calzetti et al. (2000) use a factor of 1.75 to con- 
vert from FIR to 1-1000 fim, the contribution to the 1-8 /xm 
being in fact of the order of a few percent (see also Dale et al. 
2001). 

O stars are linked not only to the galactic thermal FIR emis- 
sion through the heating of their placental clouds, but also to 
the galactic nonthermal radio emission through the acceleration 
of particles (which emit nonthermal synchrotron radiation) dur- 
ing their final SN explosion. This is reflected in a strong radio- 
FIR correlation which can be quantified (Helou et al. 1985) by 
a parameter, 



QFIR = log ( ] - l0g(/i.4GHz) 



(5) 



(with /fir in W m i^60pm=3.75xl0^2 Hz, and /i.4GHz in 
W m^^ Hz^^), that turns out to have a value 



(jFiR ^ 2.35 ±0.02 



(6) 



for local samples (Condon et al. 1991a; Yun et al. 2001; Bell 
2003). In principle then, the nonthermal radio emission pro- 
vides us with another sensitive tracer of the massive stellar 
population and hence of the instantaneous SFR (Condon 1992). 
For local SFGs, we adopt the calibration between 1.4 GHz lu- 
minosity and SFR, 

-^1.4 



SFR = 



(7) 



1.61 X 1028crgs-iHz-i 

derived by Schmitt et al. (2006) assuming a Salpeter stellar 
IMF with mass limits 0.1-100 Mq. Once adjusted for the same 
mass interval, this calibration produces SFRs a factor of ~2 
lower than the cahbration of Condon (1992), which was based 
on the Galactic relation between nonthermal 1 .4 GHz luminos- 
ity and SN rate; however, it is very similar to the more recent 
calibrations of Yun et al. (2001) and Bell (2003). (The above 
relations in eqs.(3)-(7) are, of course, self-consistent.) 

5. Analysis and results 

To evaluate L^^, we model the measured XPLFs as linear 
combinations of the 'universal' young and old XPLFs (see 
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Fig. 2. Published fits to measured XPLFs (dot-dashed lines) are modelled as the sum (solid line) of a ' young -XPLF' (dashed 
line; see eq.[l]) and an 'old-XPLF' (dotted line; see eq.[2]). In each case, the cumulative old-XPLF component is normalized at 
the break luminosity, Lbr=5x 10^^ erg s~^, and its profile can be fitted with the following expression: (1 -|- a:~^ °^)/(l -|- x^'^) 
with x=L/L\3T.. 



eqs.[l],[2]). The normalization ratio of the young to old differ- 
ential XPLFs can be expressed as a function of the fractional 
young-XP luminosity, -q, according to: 



^^v,o 

X 



1 - 

r2-/3v 



-L 



2-/3. 



2 - /3v,: 



+ 



2 - /3v,2 



r2-/3y 



L 



^min 

2-/3y- 



-/3v,i 



2-/3y 



(8) 



Since going back to the original XPLF data that are scattered 
in the literature is beyond the scope of the present work, in 
this paper we shall apply the young/old-XP decomposition to 
the published fits of the observed XPLFs. Only in the cases 
of NGC 628 and NGC 1569 did we choose to decompose the 
actual data. After performing random checks, we are confident 
that introducing this approximation has not biased the results 
of our modeling in any significant way. Our results are shown 
in Figs. 2 and 3. 

One further step involves correcting Lxp for incomplete- 
ness of the corresponding XPLFs. In nearby galaxies the 
XPLFs can be measured down to Lmin'^lO^^ erg s~^, while in 
more distant galaxies the XPLFs can only be measured down 
to higher limiting luminosities, hence Lxp are biased low with 
distance. We correct for this bias - at least approximately - by 
assuming that all XPLFs can be extrapolated with the same 7 
down to L,nin=10'^^ erg s^^. This leads to a new, distance-bias- 
corrected Lxp. The flat XPLFs (7<1, except for NGC 628) 
ensure that these extrapolations are quite reasonable. 
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Fig. 3. The measured XPLFs of NGC 628 and NGC 1569 (his- 
tograms), not shown in Fig.l, are modelled as the sum (soUd 
line) of a 'young-XPLF' (dashed line; see eq.[l]) and an 'old- 
XPLF' (dotted Une; see eq.[2]); their pubUshed power-law fits 
(see Table 1) are also shown, as arbitrarily shifted long-short- 
dashed lines. Similar information is displayed for NGC 4449 
and NGC 4631 to exempUfy our decomposition method on 
original XPLF data (as opposed to the corresponding published 
XPLF fits, shown in Fig.l). In all cases, the power-law fits are 
displayed over the luminosity ranges where they were origi- 
nally derived (see Table 1). NGC 4449 shows a clear case of 
incompleteness in source counts at low luminosities. AU data 
are from Colbert et al. (2004). 
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Fig. 4. The FIR-based SFR versus 2-10 keV luminosity re- 
lation, using the total luminosity {bottom) and the collective 
luminosity of young XPs (top) for the local sample of star- 
forming galaxies (filled circles: see Tables 1,2; NGC 3077 
is represented by a cross) and the more distant sample of 
starburst-ULIRGs (empty squares: see Table 3). In the latter set 
of high-SFR galaxies the emission is plausibly due to young 
XPs, hence the total 2-10 keV luminosities, L^, are used in 
both panels. The solid line in the top panel shows the relation 
in eq.(9); the dashed line in the bottom panel shows the relation 
ineq.(lO). 



In FigA-left the FTR-deduced SFRs are plotted versus 
L^^. The data clearly suggest a Unear correlation. 



XP 



(9) 



(with a ~20% statistical uncertainty on the calibration). Due to 
its normaUzation by a ratio of luminosities, the above relation 
is independent of the adopted cosmological model. 



6. Discussion 

In this section we discuss some issues of relevance to the vi- 
abiUty and use of L^-^^ as a SFR estimator. Specifically, we 
discuss properties of shape and calibration, and some underly- 
ing uncertainties, of the SFR-L^^^ relation. 

6. 1. Linearity 

Our SFR-L^^^ relation appears to be linear over about 5 
decades in luminosity and SFR. This result is consistent with 
that of Colbert et al. (2004) who, using a sample of 32 elhptical 
and spiral galaxies, found that Lxp was linear in both galaxy 
stellar mass and SFR (measured by A'-band and FIR+UV lumi- 
nosities, respectively). In our sample, however, poor statistics 
prevent a fuU assessment of the relation in the very-low-SFR 
regime. 

In an attempt to gain furher insight, we examined the be- 
haviour of the very-low-SFR galaxy, NGC 3077, in the L^^^- 
SFR plane. No pubUshed XPLF is available for this galaxy, but 
the measured spectra of its 6 detected XPs (Ott et al. 2003) 
can be used to single out young sources. The very soft sources 
SI, S5, and S6 are proposed by Ott et al. to be SNRs, and the 
hard source S3, speculated by Ott et al. to be an accreting bi- 
nary, is spectrally consistent (r~l) with being a HMXB. It was 
suggested by Ott et al. that S2 is an accreting binary or a back- 
ground AGN: in the former case, its r~1.65 slope may suggest 
a LMXB interpretation. The supersoft source S4 is argued by 
Ott et al. to be either a hydrogen-burning white dwarf or an iso- 
lated NS. If the identification of 4 young XPs is correct, then 
NGC 3077 approximately agrees with the SFR-L^^^ relation 
defined for our composite sample (see cross in Fig.2-left). 

NGC 2403 is another galaxy with very low SFR, poten- 
tially useful to investigate the low-L behavior of the SFR- 
L^'^^ relation. However, it seems to be a problematic object. 
As noted by Schlegel & Pannuti (2003) and discussed also 
by Fabbiano (2005), NGC 2403 is X-ray overluminous for its 
FIR luminosity and apparently violates Kilgard et al.'s (2002) 
correlation between XPLF-slope and FIR-based SFR. These 
apparent contradictions may be resolved assuming that SF in 
NGC 2403 turned off a few 10^ yr ago: the FIR emission in 
the star-forming clouds, powered by OB stars, would be dras- 
tically reduced by now, whereas the HMXBs would still be 
shining (Schlegel & Pannuti 2003). We suggest an alternative 
possibility to reconcile NGC 2403 's X-ray and FIR proper- 
ties. Spectral modelling of the 4 brightest XPs of NGC 2403 
(Schlegel & Pannuti 2003) suggests a LMXB nature for two of 
these (sources 1 and 28), and a BHXB nature for the other two 
(sources 20 and 21). For the latter we also suggest a low-mass 
donor interpretation, based on the galaxy's very low SFR and 
the lack of spatial correlation with hkely SF sites (see Schlegel 
& Pannuti 2003). If LMXB -related, tiie 4 brightest sources of 
NGC 2403 would no longer be unclassified and hence should 
be removed from the unclassified (as young or old) Ust of XPs. 
(Being old sources, these four sources would of course be un- 
suitable to trace the ongoing SF.) Of the remaining 37 less lu- 
minous unclassified XPs, distributed now in a truncated XPLF, 
12 are estimated to be interlopers (apparently distributed ran- 
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Table 3. Data II: starburst ULIRGs. 



Source 


z 






r(c) 


/12' 




Ad) 
• ' 60/j. 


Am 


FIR. 


tU) 

iK, 


Sfr(9) 


/l'4 




IRAS 12112+0305 


0.072 


293.0 


-13.79 


41.22 


0.08 


0.66 


8.18 


9.46 


45.60 


45.81 


297 






IRAS 14348-1447 


0.082 


334.5 


-13.72 


41.41 


0.07 


0.55 


6.82 


7.31 


45.62 


45.84 


315 






IRAS 15250+3609 


0.055 


222.9 


-13.64 


41.13 


0.16 


1.31 


7.10 


5.93 


45.26 


45.48 


136 


0.013 


29.89 


Arp220 


0.018 


72.3 


-12.74 


41.06 


0.61 


8.00 


104.0 


112.0 


45.48 


45.69 


225 


0.30 


30.27 


IRAS 17208-0014 


0.043 


173.8 


-13.34 


41.22 


0.20 


1.61 


9.53 


11.05 


45.21 


45.43 


122 


0.10 


30.56 


IRAS 20100-4156 


0.129 


531.6 


-13.65 


41.88 


0.1 


0.57 


5.44 


5.41 


45.90 


46.12 


595 


0.02 


30.83 


IRAS 22491-1808 


0.078 


317.6 


-14.19 


40.89 


0.06 


0.54 


5.54 


4.64 


45.46 


45.68 


216 







Luminosity distances (inMpc), computed from redshifts according to _Dl=(2c/HoS1o) [f^o^+(f^o— 2) (x/1 + SIqz— 1)] (Peacock 1999). 
2-10 keV fluxes (in erg cm~^s~^ and given in log form), from Iwasawa et al. (2001) for Arp 220 and from Franceschini et al. (2003) for the 
remaining sources. 

^'^^ 2-10 keV luminosities (in erg s~^), in log form. 

IRAS flux densities (in Jy), from Genzel et al. (1998) and Sanders et al. (1995) for 20100-4156 and from Sanders et al. (2003) for the 
remaining sources. 

FIR luminosities (in erg s~^), in log form. 

IR luminosities (Lirs1.65xLpir; in erg s"'^), in log form. 

Star-formation rates (from eq.3), in yr^^. 

1.4 GHz flux densities (in Jy), from White & Becker (1992) for IRAS 15250+3609 and from Condon et al. (1996) for the remaining sources. 
1.4 GHz luminosities (in erg s~^ Hz~^), in log form. 



domly with luminosity, see Schlegel & Pannuti 2003). Finally, 
25 sources are left as the galaxy's XP population out of which 
young sources have to be extracted by XPLF modelling. The 
young-XP luminosity estimated with our approach does com- 
ply with the SFR-L^^^ relation (see Fig.2-/e/j). Further de- 
tailed study of NGC 2403 and its XP population is clearly much 
needed, both for its intrinsic interest, and for investigating the 
lowest-L reaches of the SFR-L^"^^ relation. 

An analogous study of the use of young XPs (specifically: 
HMXBs) as SFR indicators led Grimm et al. (2003) to suggest 
a non-linear regime, SFR oc (luminosity)"-^ for SFR^ 4.5 Mq 
yr~^, in the SFR-X-ray-luminosity relation. This behaviour 
was attributed by Gilfanov et al. (2004b) to effects of low- 
numbers statistics in the distribution of XPs in low-L^;, low- 
SFR galaxies. However, we suggest that differences in our sam- 
ple selection and that of Grimm et al. may also play a role. 
The procedure of Grimm et al. involved a priori selection of 
galaxies with high SFR-to-stellar-mass ratios (based on dy- 
namical estimates for the stellar masses of galaxies, and SFRs 
derived from a variety of indicators) to ensure that would 
be HMXB-dominated and hence a tracer of the ongoing SFR. 
In contrast, we use all galaxies with measured source counts: 
by decomposing their XPLFs into young and old component, 
we a posteriori obtain the young-XP luminosity. 

Given this uncertain situation, the faint limit of the SFR- 
I/^^^ relation clearly needs further investigation. 

6.2. Calibration 

The calibration of our SFR-L^^p relation, (0.75±0.15)xl039 
erg s~^ per MQyr^^, is compatible with that, ~10^^ erg s~^ 
per M0yr~^, derived by Grimm et al. (2003) using a variety of 
SFR estimators. Griimn et al. used in their relation, but for 



their sample, which was selected such that would be largely 
HMXB-dominated (see section 6.1), L^^L^^ by construc- 
tion. 

Our calibration agrees also with the corresponding cal- 
ibration of Colbert et al.'s (2004) bilinear correlation be- 
tween Lxp and host-galaxy SFR and steUar mass. To see 
this we should account for the different definition of the vari- 
ables, SFR and X-ray luminosity, in Colbert et al.'s eq.(7) 
and in our eqs.(3),(4). Specifically, we convert their ipiR+uv- 
based definition of SFR into our adopted LiR -based definition, 
which gives SFRfir+uv=1-067SFRir, and transform their 
0.3-8 keV luminosities into our 2-10 keV luininosities (as- 
suming with Colbert et al. F = 1.7 PL spectra), which gives 
io. 3-8=1-529 L2-10- In terms of our variables, Colbert et al.'s 
bivariate relation can be rewritten as 

Lxp - (0.85 ± 0.13) x lO^^ M = (0.49 ± 0.21) x 10^^ SFR 

(with 2-10 keV XP luminosities in erg s~^, masses in Mq, and 
IR-derived SFR in Mq yr~^). Comparing this expression with 
eq.(9) one sees that: 

( ;) our young-XP luminosity corresponds to Colbert et al.'s XP 
luminosity minus a quantity, proportional to the galaxy's stellar 
mass, that represents the old-XP luminosity: LJ^^ = Lxp — 
(0.85 ±0.13) X 1G29M; and 

( ii) once adjusted for the same definitions of SFR and luminos- 
ity, Colbert et al.'s calibration, (0.49 ±0.21) x 10^^ erg s~^ per 
M0yr~^, is consistent with ours. 

6.3. Uncertainties 

The main uncertainties concern the precision with which our 
adopted young-XPLF and old-XPLF have been detemined. 
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Old-XPLF. The old-XPLF proposed by Kim & Fabbiano 
(2004), though the most updated and reliable available, extends 
in luminosity down to only ^SxlO''^ erg s"'^, whereas mea- 
sured XFLFs often reach down to ~10^® erg s~^. In our analy- 
sis we chose to extrapolate Kim & Fabbiano's function down to 
~10^^ erg s^^, but this assumption may not be fully realistic. 
A flattening of the old-XPLF at L< lO^'^ erg s~\ with a differ- 
ential slope of ~1, is suggested by counts of bulge LMXBs of 
some nearby spirals (Gilfanov 2004) and of XPs in NGC 5128 
(Voss & Gilfanov 2005). However, the situation concerning 
the low-i XPLFs may be rather compUcated. For example, 
in M Si's very well studied LMXB population clear differ- 
ences are seen between bulge and globular-cluster XPLFs, 
and different low-L breaks appear in imier-bulge, outer-bulge, 
and globular-cluster XPLFs (see Fabbiano 2005 and references 
therein). 

Our choice is motivated by concerns of unexplored com- 
plexities in the low-L old-XPLF, as well as by the considera- 
tion that in SFGs the details of the low-L old-XPLF may not 
be crucial, if Lxp is dominated by bright sources - i.e. when 
the cumulative XPLF index is 7<1, as is the case for the SFGs 
in Table 1 (except for NGC 628, whose measured XPLF how- 
ever extends over a range where Kim & Fabbiano's function is 
defined). 

In conclusion, use of Kim & Fabbiano's (2004) combined 
broken-PL old-XPLF probably represents a significant im- 
provement in accuracy over the use of a single-PL with (cu- 
mulative) slope ^ 1 as suggested by individual XPLFs (e.g., 
Kim & Fabbiano 2004; Colbert et al. 2004). We consider this 
assumption to be adequate for modeling XLPFs that do not ex- 
tend much lower than L~10^^ erg s~^. But accurate modeling 
of deeper XLPFs (either measured or extrapolated) wiU require 
better knowledge of the low-L old-XPLF. 

Young-XPLF. Either measured from very active global starburst 
galaxies (e.g.: NGC 4038/9, Zezas & Fabbiano 2002) or from 
spatially resolved young XP population (e.g.: M 81, Tennant 
et al. 2001; M 83: Soria & Wu 2003), young-XPLFs turn out 
as single PLs with (cumulative) slopes — 0.5±0.1 - the flatter 
slopes were measured in more intensely star-forming galaxies 
(Kilgard et al. 2002). Grimm et al. (2003) suggested the exis- 
tence of a 'universal' young-XPLF described as a single-PL, 
with cumulative slope —0.6 and normalization proportional to 
the SFR, in the luminosity interval ~4x 10''^— 10**" erg s^^. 
However, as discussed in Fabbiano (2005), a significant scat- 
ter of individual XPLF behaviors is observed. In our analysis 
we adopted Grimm et al.'s suggested regularity of a universal 
single-PL XPLF, but we chose the slightly flatter slope —0.5 
which we feel (following Kilgard et al. 2002) to better rep- 
resent homogeneous young-XP populations in high-SFR envi- 
ronments. 

A further complication is possible. If in a high-SFR envi- 
ronment the stellar IMF is top-heavy (e.g., Doane & Mathews 
1993; Rieke et al. 1993), the resulting correlation between SFR 
and stellar IMF implies a proportionally higher number of mas- 
sive stars, and hence a flatter young-XPLF, for higher SFRs. 
This would in principle challenge the concept of a 'universal' 
young-XPLF. If so, using one same young-XPLF to model the 



Table 4. Data 111: Hubble Deep Field North galaxies. 



Source 


z 


/l.4 




SFR 






134 


0.456 


210 


30.00 


62.1/ 112.0 


-15.55 


41.13 


136 


1.219 


180 


30.87 


455.5/821.3 


-15.72 


41.96 


188 


0.410 


83 


29.50 


19.6/ 35.3 


-16.24 


40.34 


194 


1.275 


60 


30.43 


167.8/ 302.5 


-15.70 


41.95 


246 


0.423 


36 


29.16 


9.1/ 16.3 


-16.12 


40.49 


278 


0.232 


160 


29.26 


11.3/ 20.4 


-15.80 


40.26 



All fluxes (2-10 keV: erg cm~^ s~ ^ ; 1.4 GHz: /xJy) are rest-frame and 
are taken from Ranalli et al. (2003). Radio/X-ray luminosities (mea- 
sured, respectively, in erg s"^Hz^^ and erg s~^), as well as X-ray 
fluxes, are given in log form. X-ray fluxes and luminosities were de- 
rived by modeling Chandra 0.5-8 keV counts with a power-law model 
(Ranalli et al. 2003). Star-formation rates, derived from the 1.4 GHz 
luminosity by means of eq.(7)/eq.(l 1), are expressed in Mq yr~^. 

XPLFs of galaxies with very different SF activities would be 
incorrect and would lead to a systematic bias in the analysis. In 
this scenario, the assumption of 'universal' young-XPLF im- 
plies that SF sites in galaxies should have very similar char- 
acteristics everywhere within individual galaxies and in differ- 
ent galaxies, the main difference between high- and low-SFR 
galaxies being the number and sizes, not the physical proper- 
ties, of such SF sites. 

Given the known uncertainties, we checked that our main 
results are not significantly altered by changing our adopted 
young-XPLF slope by ±0.1, which is believed to represent a 
reasonable uncertainty in the 'universal' young-XPLF slope. 

7. Star formation in thie nearby Universe 

To sample SF more completely in the nearby Universe, we 
extend our analysis to a sample of starburst-ULIRGs ^ Their 
very-high-SFRs (^ 100 M0yr~'^), and the apparent flatness 
of their 2-10 keV spectra (r~1.2: e.g., Franceschini et al. 
2003) reminiscent of Galactic HMXB spectra, suggest that 
in these objects Lx mainly originates from young XPs. Then 
L^^^~La;, and the current SFR of these galaxies is effectively 
traced by L^. Starburst-ULIRGs are natural caUbrators for the 
SFR-X-ray-luminosity relation. From FigA-left we see that, 
using Lx, our ULIRGs (empty squares) do lie on the extrap- 
olation of die SFR-L^^P relation defined (at lower SFRs) by 
the SFG sample. 

For a comparison with the SFR-L^^^ relation in eq.(9), 
we now consider the SFR-L^: correlation. For sample 1 this is 
reproduced by (see Fig.2-right, filled circles): 



SFR(>O.1M0) = 



L. 



3.8 X 



Mqyt- 



(10) 



(with a ~20% statistical uncertainty on the calibration). This 
expression is consistent with an analogous expression by 

' By this definition we mean ULIRGs that show no obvious X-ray 
spectral evidence of harboring a central AGN (e.g., Franceschini et al. 
2003). 
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-0.4 -0.2 0.2 41 42 43 44 45 46 
log (160/1100) log L(FIR) 

Fig. 5. The young- to old-XP 2-10 ke V luminosity ratio plotted 
versus the temperature index /eo/Zioo (left) and the FIR lumi- 
nosity (right) for the galaxies in sample 1 . The young-XP frac- 
tion is generally higher in more FlR-luminous galaxies, but the 
correlation appears to be tightest with the starburst phase (indi- 
cated by /go //loo)- it is highest in peak-starbursts and gets pro- 
gressively lower in evolved-starbursts and post-starbursts (see 
also Fig.l). The outlier in the right panel is NGC 1569 whose 
IRAS colors are suggestive of a peak-starburst phase in spite of 
its low luminosity. 



Ranalli et al. (2003), based on nearly the same data, once their 
different definition of SFR, involving Lfir instead of Ljr, is 
accounted for.) Inspecting the two relations in Fig.4, it is clear 
that sample 1 complies with both, whereas sample 2 does not. 

The mismatch between the SFG sample and the ULIRG 
sample in the SFR-I/^; plane is so considerable that no simple 
function can adequately characterize our full combined sample 
over the SFR range (-2^ log[SFR/(M0yr-i)];$ 2.5) consid- 
ered here. Since there is nothing special about the way our star- 
forming objects were selected, we suggest that the discrepancy 
is real and does not stem from a known bias. Our interpretation 
of the discrepancy is as follows. In all star-forming galaxies vir- 
tually the total LpiR traces the istantaneous SFR, whereas Lj., 
that is emitted partly by LMXBs (which represent the SFR of 
previous epochs) and partly by young sources (SNRs, HMXBs 
which trace the ongoing SFR), traces the integrated (from some 
previous epoch to the present) SFR. In those galaxies, there- 
fore, only a fraction of is related to the instantaneous SFR 
- a fraction which is large in SF-dominated galaxies (like our 
ULIRGs) but can be quite small in more quiescent spirals (e.g., 
the Galaxy) (see Fig.5) 

The simultaneous validity of both the SFR-_L5;;^^ and SFR- 
Lx relations for our local SFGs (with SFR<50 Mq yr~^) sug- 
gests that in these galaxies the SF activity has remained essen- 
tially constant over the past several 10^ yr, so the corresponding 
SFR (traced by L^) is, for most galaxies of sample 1, approx- 
imately the same multiple of the instantaneous SFR (traced by 
Ll^^). 

It would be misleading, however, to infer the ratio of the 
integrated to instantaneous SFR from the normalization ratio 
(~5) of eq.(lO) to eq.(9). We emphasize, in fact, that in sam- 
ple 1 there is a systematic discrepancy between the XP lumi- 
nosity (used to build Uj^^=riLxp) and the integrated lumi- 
nosity: <Lxp>— 0.40lo;|!)g<La;>. This discrepancy may orig- 
inate from inaccurate modeUng of the individual XP spectra 
(affecting Lxp; e.g., Schlegel & Pannuti 2003) and/or of the in- 



— q™=2.60 

— q„„=2.35 




27 26 29 30 31 
log L(1.4 GHz) 



26 28 30 32 

log L(1.4 GHz) 




29 30 31 32 33 
log L(1.4 GHz) 



29 30 31 32 33 
log L(1.4 GHz) 



Fig. 6. The FTR-radio correlation (left) and, equivalently, the 

FIR/radio parameter qpm. (left) for several samples of star- 
forming galaxies. Top: The local galaxies of sample 1 (filled 
circles, dotted Une) and the relatively faraway ULIRGs of sam- 
ple 2 (empty squares, solid line). The dashed line denotes 
the characteristic relation for local field galaxies according to 
Condon et al. (1991a), Yun et al. (2001), and Bell (2003). 
Middle: The high-redshift ULIRG sample of Stanford et al. 
(2000). All data are used; when 100/xm detections could not 
be obtained, la upper limits were used to compute Xfir- 
Bottom: Stanford et al.'s (2000) sample, with non-detections 
and sources with 2:>0.45 removed. 

tegrated galaxy spectra (affecting L^, e.g. Dahlem et al. 2000), 
or from the presence of a deeply buried AGN (e.g.: Delia Ceca 
et al. 2002; Komossa et al. 2003; Ballo et al. 2004) or of tinly 
diffuse emission (e.g., Griffith et al. 2000), or from a combi- 
nation of all these. Whatever its origin, the discrepancy be- 
tween Lx and Lxp contributes significantly to the difference 
between the normalizations of eqs.(9) and (10). As a com- 
parison, from Table 2 we derive <r;>~0.50±0.07 for sample 
1. Hence <L^^>'^0.2 <Lx>, as implied by the respective 
normalizations. (Incompleteness corrections, being relatively 
modest, do not substantially alter this result.) 

8. Star formation at high redshift 

Knowledge of the cosmological SF history is crucial to con- 
strain models of galaxy evolution. One key step forward in this 
direction is developing our ability to measure the ongoing SFR 
in galaxies at cosmological distances, in a way that is mostly 
unaffected by absorption. 

Sample 3 is a set of distant (-2~1) Hubble Deep Field North 
galaxies (HDFNGs) with available 1.4 GHz flux densities and 
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Chandra-hased 2-10 keV fluxes (see Table 4). No FIR data are 
available. In Table 4 we report (from RanalU et al. 2003) their 
(k-corrected: e.g., Bauer et al. 2002) rest-frame 2-10 keV lu- 
minosities and radio luminosity densities. The former are com- 
puted from Chandra counts in the soft (0.5-2 keV) and hard 
(2-8 keV) band, assuming a power-law (PL) model. If a more 
realistic model is adopted, e.g. a sub-keV thermal plus a hard 
PL model (see Dahlem et al. 1998), the resulting 2-10 keV lu- 
minosities would be shghtly (~10%) lower. 

The SFRs of our HDFNGs, estimated from ii.4 us- 
ing the (locally derived) conversion in eq.(7), are high 
(SFR~100-1000 Mgyr-^: see Table 4), suggesting that these 
galaxies are SF-dominated (see also Cohen 2003). This, and the 
consideration that at the epoch corresponding to z^l (~6 Gyr 
for our adopted cosmology) there had been no time for LMXBs 
to form, lead to the expectation that the of our HDFNGs are 
largely dominated by young XPs. We then assume L^-^^^L^. 
In the context of our current analysis, HDFNGs appear similar 
to ULIRGs. 

The link between a starburst's instantaneous SFR and ther- 
mal FIR emission is well established, rendering LpiR a rel- 
atively accurate estimator of the instantaneous SFR in star- 
burst galaxies, most notably in their peak phase (see section 
4). This is not the case for the non-thermal radio emission, 
whose calibration with the SN rate is not known precisely (see 
Condon 1992; Condon et al. 2002) and whose characteristic 
synchrotron loss timescale strongly depends on the magnetic 
field - which can be very different in galaxies of (e.g.) very 
different SFR. This means that, in principle, in a given sam- 
ple we should check the cross-correlation of the FIR-based and 
radio-based SFR indicators to ensure that the two sets of SFRs 
derived for that sample are mutually consistent. For a sam- 
ple for which the radio SFR indicator is not known directly, 
combining the FIR SFR indicator in eq.(3) with the observed 
FIR-radio correlation for that sample will yield a suitable radio- 
based SFR indicator. 

Our direct knowledge of the FIR-radio correlation 
for deep samples is still quite limited. Garrett (2002), 
analyzing a sample of distant (z^ 1.3) HDFNGs with 
28.5^S log[Li.4GHz/(crg s~^Hz~"'^)]^ 32 (i.e., overlapping in 
luminosity with sample 2), could reach no conclusive results 
because the FIR fluxes of those galaxies, which were not di- 
rectly accessible, had to be extrapolated from available ISO 
15/xm data assuming a starburst template: the resulting gpiR 
depended crucially on the adopted template, hence no infor- 
mation could be obtained on the actual value of qpin for 
that sample. Appleton et al. (2004), based on Spitzer 70/im 
data for a distant galaxy sample in the luminosity range 
10^^^ -Z^20cm/(W Hz"^)^ 10'^° (i.e., overlapping in luminos- 
ity with sample 1), concluded that g7o~2.15 out to z~2. 

Attempting to circumvent any lack of direct knowledge, 
and pushing further our assumption of similarity between 
HDFNGs and ULIRGs, we assume that any FIR-radio corre- 
lation observed for the latter wiU also be representative of the 
former. Our own data suggest an ULIRG value of gpiR— 2.6 



(Fig.6-fop). ^ The smallness of our ULIRG sample is somewhat 
compensated for by its homogeneity, as it includes only peak- 
phase starbursts (i.e., objects with /eo/Zioo"^!)- This possibly 
is the reason for the neghgible scatter of our ULIRGs around 
9fir=2.6 in Fig.6-top. A likely confirmation comes from an 
analysis of Stanford et al.'s (2000) sample of high-z ULIRGs 
(with 29.5^ log[Li.4GHz/(erg s~^Hz^^)]^ 31, i.e. spanning 
the same luminosity range as our ULIRGs): although we find 
2~ Qfir^ 2.6, in agreement with Condon et al.'s (1991b) ear- 
her result for a flux-limited IRAS sample of starbursts and 
ULIRGs, nevertheless, when obvious outUers and IRAS non- 
detections are removed, we find that most of Stanford et al.'s 
(2000) data are consistent with (7fir~2.5 — 2.6, definitely 
higher than the 'canonical' local value of -^2.35 (see Figs. 6- 
middle, 6-bottom). 

Based on the above considerations, we tentatively propose 
9fir=2.6 for the 'pure starbursts' represented by the ULIRGs 
in sample 2. If, according to our assumption, this value is 
also representative of our HDFNGs, then by combining it with 
eqs.(3)-(5) we obtain a consistent radio SFR indicator, 

oT^rj _ ^1-4 /lis, 

8.93 X 1027ergs-iHz-i ' ^ ^ 

which we suggest may apply to ideal starbursts. The proposed 
conversion happens to be close to that of Condon (1992), who 
assumed a Galactic calibration of the non-thermal radio lumi- 
nosity with SN rate, and to be higher by a factor of 1 .8 than that 
of Schmitt et al.'s (2006). We shaU use eq.(l 1) alongside eq.(7) 
to estimate SFRs from 1.4 GHz luminosities for our sample 3 
galaxies. 

Concerning the 'nonthermal-luminosity versus SN-rate' 
calibration issue, we point out that Yun & Carilli (2002), us- 
ing a FIR-radio spectral template in which the normalization 
of the non-thermal radio continuum had been left free to vary 
in order to determine a normalization most suitable for star- 
burst galaxies, found that a local sample of FIR-luminous 
(>1O^^L0) galaxies were best fit by the Galactic normaUza- 
tion adopted by Condon (1992); and that the spectral tem- 
plate, incorporating such Galactic normalization, could fit the 
observed SEDs of some distant (z^l), intensely star-forming 
(SFR-200- IOOOM0 yr-i) galaxies. So Yun & CarilH's 
(2002) results may provide circumstantial evidence that in 
very-high-SFR evironments the calibration of the radio-SFR 
conversion may be quite similar to that adopted by Condon 
(1992) - and hence in implicit agreement with our proposed 
conversion in eq.(ll). 

Using the data in Table 4, we plot in Fig. 7 the HDFNGs 
on the SFR versus X-ray luminosity plane (large filled trian- 
gles), for both L^-^^ (top) and (bottom). The left panels 
use Schmitt et al.'s (2006) conversion in eq.(7), whereas the 
middle panels use the our proposed empirical conversion in 
eq.(ll). The right panels use eq.(ll) and a two-component 
(soft-thermal plus hard-PL) model, which are probably more 
realistic for starburst galaxies than Ranalli et al.'s (2003) sim- 
ple PL model upon which the values of -Lx in Table 4 are based. 

^ Also from Fig.6, we point out that sample 1 has gFiR~2.2, hence 
it offers an essentially unbiased representation of the local SPG popu- 
lation. 
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Fig. 7. The SFR versus 2-10 ke V luminosity relations, using the collective luminosity of young XPs (top) and the total luminosity 
(bottom), for the distant sample of HDFNGs in sample 3 (large filled triangles) and, as a comparison, for the local SFGs in 
sample 1 (small filled circles) and the nearby ULIRGs in sample 2 (small empty squares). Because the emission of the HDFNGs 
is arguably due to young XPs, their total 2-10 keV luminosities, Lx, are used in both panels. These are computed from Chandra 
counts assuming a single-PL model (left, middle) and, alternatively, a soft- thermal plus hard-PL model (right). In the latter 
case, in order to maximize the effect we assumed the thermal component to have 0.2 Zq (e.g., Tsuru et al. 1997) and the PL 
component to have r=1.2 (e.g., Franceschini et al. 2003). (All new luminosities were k-corrected accordingly.) HDFNG SFRs 
were computed from 1.4 GHz luminosity densities using Schmitt et al.'s (2006) locally-calibrated conversion in eq.(7) (left) and 
our proposed 'ideal starburst' conversion in eq.(ll) (middle, right). Shown are also the relations in eq.(9) (top, soUd line) and in 
eq.(lO) (bottom, dashed line). 



Inspection of Fig. 7 leads us to conclude that: (a) if Schmitt et 
al.'s (2006) local calibration holds also at high redshifts (or, 
alternatively, for strong starbursts), then our distant HDFNGs 
don't really seem to strictly follow either relationship, falling 
somewhere in between the two; (b) if our HDFNGs comply 
with eq.(ll), then they match the SFR-L^-'^^ relation but not 
the SFR-La; relation; (c) the previous point is further strength- 
ened if the HDFNG 0.5-10 keV spectra are similar to those of 
local starburst galaxies (e.g., Dahlem et al. 1998). 



Therefore, the issue of where HDFNGs are located in the 
X-ray versus SFR plane is still unsettled. More investiga- 
tions are needed to effectively measure the ongoing SFR in 
these galaxies. Broad-band IR photometry would clearly prove 
crucial given the estabUshed role and effectiveness of the 8- 
1000 /xm luminosity as a SFR indicator. Once calibrated - us- 
ing IR-derived SFRs - on distant HDFNGs, the X-ray-based 
SFR indicator could be used to gauge ongoing SFRs in yet 
more distant galaxy samples. 



9. X-ray versus radio SFR indicators 

Finally, we plot all the galaxies of our combined samples on 
the X-ray versus radio luminosity plane (Fig. 8): the correlation 
involving Lx (left) is strong, whereas the one involving 
(right) is weak - a piecewise behavior is discernible, with a 
break at L^^^^SxlO^^ erg s~^. We interpret the different be- 
havior as follows. 

(a) The radio luminosity is a function both of the 

SN rate (hence the massive SFR) and of the aver- 
age magnetic field strength. The syncrotron loss timescale 

is ^s=(|;S^7e^^)-^^1.25xl0iO(i3;^ E^l^ yr (where 
B^Q is the magnetic field strength measured in /iG, and -EceV 
is the elecron kinetic energy measured in GeV). Assuming 
£J~5 GeV (typical of electrons radiating at ~1.4 GHz) 
and S^G~l-5 gives Ts~10*-2.5x10^ yr. Synchrotron loss 
timescales are then typically (much) longer than the typical 
SF timescale, rsF~10^ yr. This means that the current radio 
emission of a local SFG may trace the SFR integrated over the 
last ^ 10^ yr. As remarked in section 7, a similar consideration 
holds for the total 2-10 keV emission, Lx. 
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log (total) log (young XPs) 

Fig. 8. 1.4 GHz luminosity densities are plotted versus 2-10 
keV luminosities, using total luminosities (left) and collective 
young-XP luminosities (right), for the combined samples of lo- 
cal SFGs (filled circles), nearby ULIRGs (empty squares), and 
distant HDFNGs (filled triangles). The hnes are solely meant 
to guide the eye. 



(b) During a strong starburst phase we expect the prompt radio 
emission to be dramatically enhanced due to the more intense 
particle accleration (and perhaps also increased magnetic field, 
e.g. Hirashita & Hunt 2006), with a corresponding decrease in 
the synchrotron loss time. Clearly, the FIR energy density is 
much higher then (perhaps as high as Ufir simlO"^ erg cm~^ 
in the starburst region; Condon et al. 1991b). Relativistic elec- 
tron Compton energy losses are then important (with very short 



loss timescales, ric=(| 7e lo-Cgcm-^ ) ^^x^"' y"' 
for the same parameters as in the previous paragraph). [Of 
course, Compton losses of synchro-emitting electrons off CMB 
photons increase with redshift as (1 -|- z)"^.] These three factors 
(enhanced fresh radio emission and shorter synchrotron- and 
Compton-loss timescales) concur in making the observed ra- 
dio emission of strong starbursts an accurate measure of their 
instantaneous SFR. 

The varying ability of synchrotron radio emission to trace 
the instantaneous SFR, as a function of the different physical 
conditions prevailing in different galaxies, may play some role 
in the discussion on whether the value of qfir is 'universal' and 
on the real meaning and caUbration of the radio SFR indicator 
(see section 8). 

In either case, the galactic radio and (total) X-ray luminosi- 
ties should measure, at any given phase of a galaxy's SF history, 
essentially the same (definition of) SFR (i.e., instantaneous in 
strong starbursts, and integrated over the past ^ 10^ yr in more 
quiescent disks). We thus expect a Unear Liaghz-Lx correla- 
tion. 

The integrated past SFR may in principle be unrelated to 
the current SFR, traced by . In our local SFGs galax- 
ies, however, if the SF activity has remained essentially con- 
stant over the past several ~10* yr (as argued above), the 
past integral SFR traces (to a multiphcative factor) the cur- 
rent SFR. This makes the integrated vs instantaneous SFR cor- 
relation possible, in the form of a linear correlation in the 
Li.4GHz,L^^ plane (Fig.S-right, filled circles). For strong 
starbursts a Unear correlation exists, too, but it is displaced from 
the SFG one given the different definition of the SFR being 
measured by the radio emission in the two cases (instantaneous 
for starbursts, integrated for SFGs). Therefore the linear corre- 



lation between young-XP emission and radio emission would 
exist separately for mildly star-forming galaxies and for ex- 
treme starbursts. Its piecewise linearity, as well as the meaning 
of measuring the current versus integrated SFR, makes L^^- 
L\A GHz a companion relation to SFR-ix- 

10. Summary and conclusion 

In this paper we examined the issue of the collective 2-10 keV 
emission of young XPs and its role as a SFR indicator in star- 
forming galaxies. For a sample of local star-forming galaxies 
with available data, we estimated the young-XP luminosity by 
modelling the observed XPLFs in terms of 'universal' young- 
and old-XPLFs. 

For SFR^ 50 MQyr"^ galaxies, the collective emission of 
young XPs turns out to correlate linearly with the FTR-based 
SFR. In this SFR range both the SFR-L^^p relation, pre- 
sented here, and the SFR-La^ relation, proposed by Ranalli et 
al. (2003), are valid. 

The relation is extended to higher SFRs by using a sam- 
ple of extremely starburst-dominated ULIRGs. Since their L^. 
is arguably dominated by young XPs and hence by the in- 
stantaneous SFR, when Lrf^^=Lx these ULIRGs comply with 
the SFR-LJ^^ relation. Similar considerations may hold for 
a sample of intensely star-forming Hubble Deep Field 
North galaxies - especially so if their radio SFR indicator is 
calibrated shghtly higher than for local galaxies, as suggested 
by the FIR-radio relation for starburst-dominated ULIRGs. 

Overall, the SFR-L^^^ relation is suggested to hold over 
the broad SFR range 0.01;$ SFR/(MQyr-i)< 1000, according 
to: 

ryXP 

^^^O = (0.75±0.15)^xl03^ergs-^ " 

For galaxies of very different SFRs, the same SFR-X-ray- 
relation holds only when the young-XP emission, not the total 
luminosity, is used. This is so because the FIR emission and the 
young-XP emission both trace the instantaneous SFR, whereas 
the total X-ray luminosity traces the instantaneous SFR in ex- 
treme starburst galaxies and the integrated SFR in more quies- 
cent ones. Consequently, a SFR-ix relation exists separately 
for very-low and very-high SFR galaxies (e.g., our local SFGs 
and distant ULIRGs, respectively). In particular, the existence 
of the SFR-_Lx relation for our local, low-SFR galaxies sug- 
gests that SF has not changed dramatically over the last ~10^ 
yr. 

The SFR-_L^^^ relation represents the most adequate X- 
ray estimator of the instantaneous SFR in galaxies. Its reflects 
the equivalence of two complementary measures of the cur- 
rent SFR, based on the observed manifestations of massive 
stars at their birth (short-lived FIR emission from placental dust 
clouds) and near their death as compact renmants (short-Uved 
X-ray emission from close binary accretors), respectively. 
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